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Summary. — Stochastic background of gravitational waves (GW) generated by 
the interactions between primordial black holes (PBH) in the early universe and 
by PBH evaporation is considered. If PBHs dominated in the cosmological energy 
density prior to their evaporation, GWs from the earlier stages (e.g. inflation) 
would be noticeably diluted. On the other hand, at the PBH dominance period 
they could form dense clusters where PBH binary formation might be significant. 
These binaries would be efficient sources of the gravitational waves. 
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The registration of the gravitational waves generated in the early universe could bring 
an important information about inflation, possible (first order) cosmological phase transi- 
tions, topological defects, such as cosmic strings, etc. Very sensitive GW detectors such 
as LIGO and LISA may make make the discovery opening a new era of gravitational 
wave astronomy. In particular, an observation of stochastic cosmological background of 
low frequency GWs could be a final proof of inflation. However, an absence of a such 
background would not mean that the universe was not in inflationary stage. First, the 
present day density of GWs depends upon the model of inflation and, second, there could 
be a mechanism which suppresses the density of GWs at post-inflationary stage. Such a 
mechanism is described in my talk. Though the inflationary background of GWs could 
be noticeably suppressed, a new higher frequency GWs would be generated by the sug- 
gested mechanism. The talk is based on two papers [1, 2]. In the second one a detailed 
reference list is presented which is reduced here due to lack of space. 

We consider GWs produced by the interactions between primordial black holes (PBH) , 
as well as by their evaporation. PBHs are supposed to be very light, so they evaporated 
before the big bang nucleosynthesis (BBN) leaving no trace in the present day universe, 
except for GWs. The life-time of an evaporating black hole with mass M is equal to [3] 
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where mpi = 1.22 • 10 19 GeV = 2.176 -10 5 g is the Planck mass and N e ff is the number 
of particle species with masses smaller than the black hole temperature: 



(2) T BH 



8ttM 



The corrections due to the propagation and back-capture of the evaporated particles [4] , 
the so called grey factor, change this result by a factor of order unity and are not included 
here. 

According to ref. [5], to avoid a conflict with BBN the life-time of PBHs should be 
shorter than t w 10 -2 s and thus the black holes should be lighter than 



(3) M< 1.75 -10 



100 



The temperature of such PBHs exceeds 3 • 10 4 GeV and correspondingly N e ff > 10 2 . 

Formation of PBHs from primordial density perturbations in the early Universe was 
considered in pioneering papers [6, 7]. PBHs were formed when the density contrast, 
Sp/p, at horizon was of the order unity or, in other words, when the Schwarzschild 
radius of the perturbation was of the order of the horizon scale. If PBHs were created 
at the radiation dominated stage, when the cosmological energy density was p(t) = 
3m 2 pl /(32Trt 2 ), and the hor izon was lh = 2i, the mass of such PBHs would be: 

(4) M(t) = m 2 Pl t ~ 4 ■ 10 38 

where t is the cosmological time. 

The fraction, fi p , of the cosmological energy density of PBH produced by this mech- 
anism depends upon the spectrum of the primordial density perturbations. If the usual 
flat Harrison- Zeldovich spectrum is assumed, then tt p would be quite small. We have not 
calculated f2 p but have taken it as a free parameter of the model. One reason for that is 
that the spectrum of the density perturbations at small wave lengths is unknown. More- 
over, there could be other mechanisms of PBH formation. In particular, in refs. [8, 9] 
a model of PBH formation has been proposed which might lead to considerably larger 
probability of PBH formation. The mass spectrum of PBHs produced by the latter 
mechanism has the log-normal form: 



dN n 

(5) dM= Cex P 



(M - M ) 2 



where C, Mq, and Mi are some model dependent parameters. Quite naturally the central 
value of PBH mass distribution may be in the desired range Mo < 10 9 g. In this model 
the value of Q p may be much larger than in the conventional model based on the flat 
spectrum of the primordial fluctuations. We will not further speculate on the value of 
fl p and on the form of the mass spectrum of PBH. In what follows we assume for an 
order of magnitude estimate that the spectrum is well localized near some fixed mass 
M and that O p is an arbitrary parameter. Different mechanisms of PBH production are 
reviewed e.g. in ref. [10]. 
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The energy density of nonrelativistic PBHs drops down as 1/a 3 , while the energy 
density of the initially dominant relativistic matter drops as 1/a 4 , where a = a(t) is 
the cosmological scale factor. So the relative contribution of PBH into the total energy 
density rises as 

(6) n BH (t) = n p a(t)/a p = n p (t/t p )^ 2 , 

here t p = M/m 2 Pl is the PBH production time, related to their mass by cq. (4). So ini- 
tially flBH{t) rises as t 1 / 2 and at some stage it reaches unity and after that flsH remains 
constant till the PBH evaporation. PBHs would begin to dominate in the cosmological 
energy density at t = t eq = M(m Pl Clp), if t eq > tbh, cq. (1). This can be transformed 
into the lower limit on the PBH mass: 

N eff\ 1/2 mpi 



(7) M>6.10-(^j 



If condition (7) was fulfilled, the universe expansion regime was initially relativistic, 
radiation dominated (RD), then after t = t eq it became non-relativistic, matter domi- 
nated (MD). Later after PBH evaporation, t > tbh the universe returned to RD stage 
again, and only after very long time, t = t^ss ~ 10 5 years, the expansion became mat- 
ter dominated. After that time the large scale structures (galaxies, their clusters, etc) 
began to form. As is known, cosmological structure formation took place at MD stage, 
when initially small primordial density perturbations started to rise due to gravitational 
instability. In our case the density perturbations started to rise at t > t eq . According to 
the theory, at MD stage A = 5p/p ~ a(t), till the perturbations remains small, A <C 1. 
When A reaches unity, the perturbations quickly rise and as a result they becomes quite 
large, A>1. In the present day universe A ~ 10 5 at the galactic scale. 

In our scenario we expect formation of high density clusters of PBHs with density 
contrast which rose as A(t) = A,„(t/tj rl ) 2 / 3 , where tt n > t eq is the moment when the 
perturbation comes inside the cosmological horizon. The density contrast would reach 
unity at t\(ti n ) such that: 

(8) A[ti (tjn)] = A m [i!(t m )/i m ] 2/3 = 1 or ti(t in ) = t in Ar 3/2 . 

To this end the PBH life-time should be longer than t\. 

After the density contrast has reached unity, the cluster would decouple from the 
general cosmological expansion. In other words, the cluster stopped to expand together 
with the universe and, on the opposite, it would begun to shrink when gravity took over 
the free streaming of PBHs. So the cluster size would drop down and both ubh and 
Pb would rise. The density contrast would quickly rise from unity to A& = Pb/p c 3> 1, 
where p c and pb are respectively the average cosmological energy density and the density 
of PBHs in the cluster (bunch). It looks reasonable that the density contrast of the 
evolved cluster could rise up to A& = 10 5 — 10 6 , as in the contemporary galaxies. After 
the size of the cluster stabilized, the number density of PBH, ubh, as well as their mass 
density, pbh, would be constant too. But the density contrast, A& would continue to 
rise as (t/ti) 2 because p c drops down as l/t 2 . From time t = t\ to t = tbh the density 
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contrast would additionally rise by the factor: 



(9) 



A{t bh ) = A b 



( TBH 



This rise is associated with the drop of the average cosmological energy density, p ~ 1 /t 2 , 
but not with the absolute rise of 6 p. This effect is absent in the present day universe 
because the time when A reached unity was close to the present universe age. 

GWs could be generated in the processes of PBH scattering in the high density clusters 
and, in particular, the GW emission could proceed from the PBH binaries. Both processes 
are strongly enhanced in the clusters. The probability of scattering and binary formation 
rate are proportional to the square of the number density of PBHs, ubh- However, the 
net effect on the cosmological energy density of the emitted GWs is linear in ubh because 
it is normalized to the total cosmological energy density. 

The cross-section of the graviton bremsstrahlung was calculated in ref. [11] for the 
case of two spineless particles (here black holes) with masses m and M under assumption 
that m -C -A/. In non-rclativistic approximation, the differential cross section is: 



(10) 



dcr 



64M 2 m 2 d£ 
15^5 T 



where £ is the ratio of the emitted graviton frequency, uj = 27r/, to the kinetic energy of 
the incident black hole, i.e. £ = 2mw/p 2 . In what follows we will use this expression for 
a simple estimate assuming that it is approximately valid for m ~ M . 

The energy density of gravitational waves emitted at the time interval t and t + dt in 
the frequency range u and uj + duj is given by 



(11) 



— = Vrel n BH ui- dt, 



where ubh is the number density of PBH and v re i is their relative velocity. The latter is 
close to the virial velocity of PBHs in the cluster and can be about 0.1. As noted by the 
authors of ref. [11], Weizsacker- Williams approximation is not valid. This means that 
there could be some difference between classical and quantum graviton emission. 

The graviton bremsstrahlung proceeded till the PBH evaporation. Hence to find the 
total energy of the produced gravitons we need to integrate their energy spectrum over 
frequencies and redshift from tbh down to the moment of the cluster formation. Thus 
we obtain for the cosmological energy fraction of GWs: 



(12) 



(breri 
GW 



{WmaxjTBH) 



N, 



eff 



100 



M ) 



Here coefficient Q > 1 reflects the uncertainty in the cross-section due to the unaccounted 
for Sommerfeld enhancement [12]. Note that A may be considerably larger than 10 5 . 
For the details see ref. [2]. 

The frequency /» of GW produced at time £» during PBH evaporation, is redshifted 
down to the present day value, f, as: 



(13) 



/ = /* 



a(i*) 



ao 
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where To = 2.725 K is the temperature of the cosmic microwave background radiation 
at the present time, T* = T(t») is the plasma temperature at the moment of radiation of 
the gravitational waves, and gs(T*) is the number of species contributing to the entropy 
of the primeval plasma at temperature T* . It is convenient to express To in frequency 
units, T = 2.7 K = 5.4 ■ 10 10 Hz. 

The density parameter of the gravitational waves at the present time is related to 
cosmological time t* as: 



(14) fW(io) = V GW (U) 



a (t*) \ ( H t 



a(t )J \H 



where H = lOO/io km/s/Mpc is the Hubble parameter and ho = 0.74 ± 0.04 [13]. 

Using expression for redshift (13) and taking the emission time f* = tbh wc obtain: 

/ inn \ 1/3 

(15) n GW (t ) = 1.67 x 10- 5 fe 2 jrp \ ) ^gw(tbh). 

Now we find that the total density parameter of gravitational waves integrated up to the 
maximum frequency is: 

10 5 g 



M 



(16) h 2 n GW (t ) «0.6-10~ 21 /v 

where K is the numerical coefficient: 

A ^ ( N eff \ ( Q \ ( 100 



(17) jr=p2f)(* 
y ' Ko.iJ V io 5 



100 J VlOOy \gs{T B H 



1/3 



Presumably K is of order unity but since A may be much larger than 10 5 , see eq. (9), 
K may also be large. 

Classical emission of GW at the scattering of non-relativistic bodies is well described 
in quadrupole approximation. If the minimal distance between the bodies is larger than 
their gravitational radii, the energy of gravitational waves emitted in a single scattering 
process is equal to: 

, 10 , rrn 37tt M 2 m 2 v 

(18) «^ = ___£_, v «l, 

where b is the impact parameter. 

The differential cross-section of the gravitational scattering of two PBHs in non- 
relativistic regime, q 2 <C 2M 2 , is: 

M 2 dq 2 2M 2 1 „ 
19) dcr = -V = -^-bdb 

m pi r m pi 
and the rate of the energy emission by the GWs is given by: 

(20) dpGHr= __p M _d*. 
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The energy density parameter of GW at the moment of BH evaporation can be obtained 
integrating this expression over time and frequency. Thus we obtain: 

(21) «w(^,-,. W -(^(^)(^)(^). 

If we allow for b ~ r g , the energy density of GWs at the moment of PBHs evaporation 
might be comparable to unity. 

Now we can calculate the relative energy density of GWs per logarithmic frequency 
interval at the present time: 

(22) ^.I^^.^^)^, 



where a' is the coefficient at least of order of unity: 




It may be much larger, if ^> 10 5 . 

More efficient mechanism of GW emission may be radiation from the PBH binaries, if 
their number in the high density clusters is sufficiently high. To form the binary bound 
state PBHs should sufficiently cool down losing their kinetic energy. The cooling could be 
achieved by the energy loss to the gravitational wave radiation discussed above and by the 
dynamical friction [14]. A particle moving in the cloud of other particles would transfer 
its energy to these particles due to their gravitational interaction. However, one should 
keep in mind that the case of dynamical friction is essentially different from the energy 
loss due to gravitational radiation. In the latter case the energy leaks out of the system 
cooling it down, while dynamical friction does not change the total energy of the cluster. 
Nevertheless a particular pair of black holes moving toward each other with acceleration 
may transmit their energy to the rest of the system and became gravitationally captured 
forming a binary. 

The dynamical friction time was estimated in ref. [15]. In both cases v > a and v < a, 
where a is the velocity dispersion, the characteristic time was of the order of: 

/ 100 \ f M\ /10 6 \ 

For PBH masses below a few grams dynamical friction would be an efficient mechanism 
of PBH cooling leading to frequent binary formation. Moreover, dynamical friction could 
result in the collapse of small PBHs into much heavier black hole. Even the whole high 
density cluster of PBHs could form a single black hole. These processes of heavier black 
hole formation would be accompanied by a strong burst of gravitational radiation. 

The emission of GWs from a binary results in the energy loss which is compensated by 
a decrease of the radius of the binary and of the rotation period. As a result the system 
goes into the so called inspiral regime. Ultimately the two rotating bodies coalesce and 
produce a burst of the gravitational waves. To reach this stage the characteristic time 
of the coalescence should be shorter than the life-time of the system. In our case it is 



(24) 



tdf 
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the life-time of PBH with respect to the evaporation. The coalescence time of the binary 
made of two BH with masses Mi and M 2 can be easily calculated, see e.g. book [16]: 



(25) 



5i?Q m 



Pi 



256M 1 M 2 (M 1 + M 2 ) 



where Rq is the initial radius of the binary. This result is true for a circular orbit of the 
binary. In the case of elliptic orbit the eccentricity drops down due to GW emission and 
the system approaches to the circular one. We may use eq. (25) for an order of estimate 
of the life-time of the binary. 

There are two interesting limiting cases, when t co 3> tbh and vice versa. In the first 
case the stationary orbit approximation is valid and each binary emits GWs with fixed 
frequency equal to twice the orbital frequency and the frequency spectrum is determined 
by the distribution of the binaries on their radius. As is shown in rcf. [2], if the stationary 
regime was realized, the spectral density parameter today would be: 



(26) 



10- 8 e 



\Neff~ 
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100 





100 
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where e is the fraction of binaries with respect to the total number of PBHs in the cluster 
and gs(T(r bh) is the number of the entropy degrees of freedom at the moment of PBH 
evaporation when the plasma temperature was equal to T(tbh)- Here the possibly weak 
redshift dilution of GWs by the factor (t co /tbh) 2 ^ 9 is neglected. 

If the system goes to the inspiral phase, then we would expect today a continuous 
spectrum in the range from f m in ~ 10 7 Hz to f max ~ 3 • 10 14 Hz. However if we take into 
account the redshift of the early formed binaries from the moment of their formation to 
the PBH decay, the lower value of the frequency may move to about 1 Hz. 

PBHs could also directly produce gravitons by evaporation. The total energy emitted 
by BH per unit time and frequency oj (energy) of the emitted particles, is approximately 
given by the equation (see, e. g. book [17]): 



(27) 



(— ) 

\dtdu J 



2N, 



TT ™ 4 



m Pi c^/Tbh - 1 



where T is the BH temperature (2). Due to the impact of the gravitational field of BH 
on the propagation of the evaporated particles, their spectrum is distorted [4] by the so 
called grey factor g(uS), but we disregard it in what follows. 

The frequency spectrum of the evaporated gravitons is not thermal because of the 
different redshifts in the course of the evaporation. According to the calculations of 
ref. [2] the spectral density parameter of GWs at t = tbh is equal to: 



(28) 

where 

(29) 



^Giv( w *; T BH ) 



\TbhJ 



2.9 • lO'WV 



"Pi 
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dz(l + z) 1/2 



exp [{z + 1)u;*/Tbh] - 1 
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and 



With respect to the thermal spectrum, spectrum (28) has more power at small frequencies 
due to redshift of higher frequencies into lower band and less power at high cj* . 
The spectral density of the evaporated gravitons today would be: 



(31) SW(/;*o)=2.7.10--^ [-jf ) I y 
where Tq is the BH temperature redshiftcd to the present time: 



(32) T = 4.5 ■ 10 ib Hz 



100 \ 1/12 / 100 \ 1/2 / M ^ 1/2 



9s{T B h)J \NeffJ \W 5 g 



The mechanisms of GWs generation considered here could create quite high cosmolog- 
ical fraction of the energy density of the relic gravitational waves at very high frequencies. 
Unfortunately at the lower part of the spectrum VLqw significantly drops down making 
such GWs outside the reach of LISA or LIGO. Still the planned interferometers DE- 
CIGO/BBO and detectors based on the resonance graviton-photon transformation could 
be sensitive to the predicted high frequency GWs. 
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